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We study the fidelity approach to quantum phase transitions (QPTs) and apply it to general 
thermal phase transitions (PTs). We analyze two particular cases: the Stoner-Hubbard itinerant 
electron model of magnetism and the BCS theory of superconductivity. In both cases we show that 
the sudden drop of the mixed state fidelity marks the line of the phase transition. We conduct a 
detailed analysis of the general case of systems given by mutually commuting Hamiltonians, where 
the non-analyticity of the fidelity is directly related to the non-analyticity of the relevant response 
functions (susceptibility and heat capacity), for the case of symmetry-breaking transitions. Further, 
on the case of BCS theory of superconductivity, given by mutually non-commuting Hamiltonians, we 
analyze the structure of the system's eigenvectors in the vicinity of the line of the phase transition 
showing that their sudden change is quantified by the emergence of a generically non-trivial Uhlmann 
mixed state geometric phase. 
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One of the main characteristics of quantum mechanics that makes it different from any classical physical theory 
is that in quantum mechanics two different quantum states, being either pure or mixed, are in general not fully 
distinguishable. By fully distinguishable we mean that it is possible, upon a result of a single-shot measurement of a 
suitable observable, to infer with probability one in which of the two given quantum states the observed system has been 
prepared. In particular, two pure quantum states are fully distinguishable if and only if they are orthogonal to each 
other. Otherwise, the maximal probability to unambiguously distinguish between two non-orthogonal pure quantum 
states is always strictly smaller than one. The reason for this lies in the fact that, while the outcomes of measurements 
on classical systems are, at least in principle, given with certainty, quantum measurements in general generate non- 
trivial probability distributions. This feature of quantum mechanics has found numerous applications within the 
field of quantum information and computation, in particular in quantum cryptography, quantum communication 
complexity, designing novel quantum algorithms, etc. (for an overview, see [H). 

Within the field of quantum information, the function widely used to quantify the distinguishability between two 
quantum states p\ and p2 is fidelity given by the expression: 



F(pi,p 2 ) = Try VfkfoVfa- (1) 

Note that in the case of pure states p\ — |-0i)(V'i| an d f>2 — IV^XV^I, the above expression reduces to 
F(\ipi){t()i\,\ip2){ip2\) — | (^l | ^2) 1 5 which is nothing but the square root of the probability for a system in state 
^2} to pass the test of being in state l^i)- The fidelity (fTJ) between two quantum states, given for two systems 1 and 
2, quantifies the statistical distinguishability between them, in a sense of classical statistical distinguishability between 
the probability distributions obtained by measuring an optimal observable in states p\ and p2- In other words, for 

every observable A, we have F{p u p 2 ) < F c {{ Pl {i\A)} , {p 2 (i\A)}) = J2i \J Pi{i\A)p 2 {i\A) , where {p a (i\A)}, a <E {1, 2}, 

is a probability distribution obtained measuring the observable A in the state p a , and F c ({pi(i\A)} , {p2(i\A)}) is the 
classical fidelity between the two probability distributions {pi(i|A)} and {p 2 (i|^4)}- For an overview of the results on 
distinguishability between quantum states and its applications to the field of quantum information, see Q and the 
references therein. 

Quantum mechanics was originally developed to describe the behavior of microscopic systems. Therefore, most of 
its applications are focused on the study of the properties and dynamics of quantum states referring to such systems, 
where quantum features dominate. On the other hand, it is a common assumption that classical behavior emerges in 
the thermodynamical limit, when the number of degrees of freedom of the system becomes large (in the limit of a large 
number of microscopic sub-systems). Yet, there is no special objection why quantum mechanics should not be generally 
applicable, even to macroscopic systems. Indeed, macroscopic phenomena such as magnetism, superconductivity or 
superfluidity, to name just a few, can only be explained by using the rules of quantum mechanics. As in these, as well 
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as in many other, highly physical relevant cases, the macroscopic features of matter are given through the features 
of its quantum states, the question of quantifying those macroscopic properties given by many-body quantum states 
arises as a relevant problem in physics. 

In this study, we are interested in those macroscopic features of matter that define its thermodynamical phase. 
Different phases of matter are separated by the so-called regions of criticality, regions in parameter space where the 
system's free energy becomes non-analytic. As the free energy is a function of the system's states, it is precisely the 
features of its states that determine the phase the system is in. Indeed, different phases have different values of the 
order parameter, given by the expectation value of a certain observable. As in the case of general quantum states, 
here as well, fidelity can be used as a function whose behavior can mark the regions of criticality (and therefore the 
phase transitions). 

In the case of quantum phase transitions (QPTs) [4J, which occur at zero temperature and are driven by purely 
quantum fluctuations, the study of the ground state fidelity has been first conducted on the examples of the Dicke 
and XY models Note that in this case, the ground states are pure quantum states, so the fidelity is given by a 
simple overlap between two pure states. It was shown that approaching the regions of criticality the fidelity between 
two neighboring ground states exhibits a dramatic drop . Subsequently, the fidelity approach to QPTs was applied 
to free Fermi systems and graphs @ and to the Bose-Hubbard model 0]. The connection between fidelity, scaling 
behavior in QPTs and the renormalization group flows was introduced in Q and further discussed in Q. Also, it 
was shown that the fidelity can mark the regions of criticality in systems whose QPTs can not be described in terms 
of Landau-Ginzburg- Wilson (LGW) theory: when the ground states are given by the matrix product states [loj ] . 
in the case of a topologically ordered QPT [llj], and in the Kosterlitz-Thouless type of transition [13] as well. The 
formal differential-geometry description of the fidelity approach to QPTs was first introduced in [l3| and subsequently 
developed in 14 1, where the connection to the Berry phase approach to QPTs (see [lj|) was also established. Further, 
on the example of the spin one-half XX Z Heisenberg chain, it was shown that the fidelity does not necessarily exhibit 
a dramatic drop at the critical point, but that the proper finite-size scaling analysis allows for correct identification 
of the QPT, see [3]. An interesting example of a Heisenberg chain where the fidelity approach fails when applied to 
the ground states, but does mark the point of criticality when applied to the first excited states, was discussed in [lil ] 
(note that the numerical results were obtained for up to 12 spins only, which leaves open the question of the ground 
state fidelity behaviour in the thermodynamic limit). Introducing the temperature as an additional parameter, QPTs 
were studied in [17j | and it was shown that extending the fidelity approach to general mixed (thermal) states can still 
mark the regions of criticality as well as the cross-over regions at finite temperatures. Finally, the genuine thermal 
PTs were discussed in [3] and [3], where the connection between the singularities in fidelity and specific heat or 
magnetic susceptibility was explicitly shown for the cases of systems given by mutually commuting Hamiltonians and 
symmetry-breaking PTs of LGW type. 

In this paper, we apply the fidelity approach to general thermal phase transitions [27| . We analyze in detail two 
particular examples given by the Stoner-Hubbard model for magnetism and the BCS theory of superconductivity. 
In general, a system is defined by a Hamiltonian H(U) which is a function of a set of parameters representing the 
interaction coupling constants generically denoted as U. In thermal equilibrium, a system's state is given by a density 
operator p(T,U). Thus, in discussing the general, thermal as well as quantum phase transitions, we can consider 
the coupling constant (s) U and the temperature T to form a "generalized" parameter q = (T,U). We consider 
the behavior of the fidelity F(p(q), p(q)) between two equilibrium thermal states p(q) and p(q) defined by two close 
parameter points q — q(T, U) and q = q + Sq = (T + 5T, U + 5U). We show that in both models considered, the PTs are 
marked by the sudden drop of fidelity in the vicinity of regions of criticality - a signature of enhanced distinguishability 
between two quantum states defining two different phases of matter, based on both short-range microscopic as well 
as long-range macroscopic features. For the general case of mutually commuting Hamiltonians, we analytically prove 
in detail that the same holds for PTs which fall within the symmetry-breaking paradigm described by the LGW 
theory (see also [Hi ] and [lil]). Further, for the case of mutually non-commuting Hamiltonians, on the example of 
BCS theory of superconductivity we show that the non-analyticity of the fidelity is accompanied by the emergence of 
a generically non-trivial Uhlmann geometric phase [l9j, th e mixed-state generalization of the Berry geometric phase 
(for the relation between QPTs and Berry phases, see [lij], [l3| and [lil]). 
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STONER-HUBBARD ITINERANT ELECTRON MODEL FOR MAGNETISM 

First, we discuss the case of the Sterner- Hubbard model for itinerant electrons on a lattice given by the Hamiltonian 

H SH = J2 £k (4 T a feT + + u^2 4i d n d li d n- ( 2 ) 

fc I 

The anti-commuting fermionic operators ct represent the free-electron momentum Bloch modes [a £ {Til}); while 
the on-site operators cj = V^ 1 / 2 J^k e ~ lkx '^k<j are gi ven by their Fourier transforms, where x\ represents the position 
of the i-th lattice site. The coupling constant U > defines the on-site electron Coulomb repulsion. There are in 
total N electrons and they occupy the volume V (such that in the thermodynamic limit, when N, V — > oo, we have 
N/V — > const). Finally, we assume for simplicity that the kinetic energy is given by Sk — fi 2 fc 2 /(2m), while the 
number operators are n kcr = ^ ka c ka and n; CT = ej^c/c = V^ 1 J2 q k e%qXl ^ka^k+qa ■ Thus, 

Hsh — y^£fc»fca- + U y^n; T n;|. (3) 

ka I 

In order to obtain the mean-field effective Hamiltonian, we neglect the term quadratic in the fluctuations 8fii a = 
fiia — ni a , where ni a — in the potential written as hi^hii — 8hi-\5hii + n^hn + n;jn;j — rijfTijj. Expressing the 

number operators in terms of the Bloch momentum operators, and using (c kcr Ck'cr) = S kk '(h ka ) (which expresses the 
translational invariance in a ferromagnetic ground state), the mean- field linearized effective Hamiltonian becomes: 

Hsh = (^T"fcT + Ekl^u) ~ VUn m . (4) 
fc 

Here, n a — N a /V is the density of electrons with spin projection along z-axis given by a € {f, 1} (N a = X)fc(^fco-) 
being the total number of electrons with spin a) . The one-particle electron energies in this effective model are obtained 
by shifting the free-electron energies Ek by an amount depending on the particle's spin: 

E k ] = £fc + Un l , 

E kl =e k + Unf. (5) 

Since the one-particle energy modes are decoupled, the overall ground state is obtained by filling the electrons up to 
the Fermi level £f'- 

\g) = ®fc<fc FT 4t ®k<k Fl 4|l°)) ( 6 ) 

where |0) represents the vacuum state with no electrons and kp-t the maximal value of the momentum for spin up 
electrons, given by Ek Fr — £f (and analogously for kpi) 28]. Note that, due to the different dispersion formulas 
([5]) for particles with spin up and spin down, in general the values of kp-\ and kp[ that minimize the ground state 
energy and therefore define the state ([6]) itself, are different and consequently the number of up and down electrons 
will be different. This is precisely the reason for the existence of magnetism in this model. As soon as the energy of 
the "biased" (magnetic) state, for which for example fc^T > kpi, becomes lower than the energy of the "balanced" 
(paramagnetic) state (fc^f = fcFj.), a magnetic phase transition will occur. Obviously, for reasons of symmetry, the 
magnetic state can be reversed with the > kp^ and fepj < kpi cases having the same energy, in the absence of 
an external symmetry breaking field H. The qualitative picture of the emergence of magnetic features can be seen 
already from looking at the original Hamiltonian @: in the U — ► limit, when the Coulomb interaction is negligible, 
the Hamiltonian represents a system of free electrons that exhibits no magnetic order (all Bloch states are doubly 
occupied); in the opposite U — > oo limit, the second term of the Hamiltonian becomes the dominant one and is 
minimized by one of two possible states for which either = or = 0. 

The quantitative analysis of the zero-temperature critical behavior of the effective Hamiltonian (j4| can be done 
by looking at the divergence of the magnetic susceptibility, using the one-particle energy dispersion relations ^ . At 
T = it leads to the well known Stoner criterion for the emergence of magnetism (20| : 



D F U C = I, 



(7) 
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where U c is the critical value of the coupling constant above which the system is in a magnetic phase and Dp = D(e F ) 
is the density of states around the Fermi energy. From this, one can obtain (see Appendix 1) the critical value U c of 
the coupling constant [29| : 

4 V 

U c = - e F — , (8) 
3 TV v ; 

where N is the total number of electrons, V is the volume of the system and m is the electron mass. 

Alternatively, the above result can be derived by minimizing the overall ground state energy, thus obtaining the 
explicit dependencies kp] — k F ^(U) and kw \ = kpi(U) which determine the ground state ([6]) and U c in particular (the 
maximum value of U for which kp-\ = kpi) [301 ] . The total number of electrons in the system is given by N = TVj + N± , 
and thus: 

^ = 6^[(*n) 3 + (*n) 3 ], (9) 

since the up and down electrons occupy spheres of radius kp-\ and kp[ , with volumes V(kp-[) and V(k F i), in momentum 
space, with a density of states V/{2ti) 3 , as follows from the periodic boundary conditions for the Bloch functions. 
From Ek F1 — Ef. F , (see Appendix 1), using the energy dispersion formulas (O, we obtain the second equation that 
determines the Fermi momenta kp-\ and kpi (with a — ^k F U c ): 

(fc FT - k Fl )[{k F] +k F i)- — (fc| T + k F] k Fl + k Fi )\ = 0. (10) 
We see that the above two equations always have the trivial solution k F ^ = k F [ = k F , which from ((9]) is k F = 

1 /3 

(37r 2 ^) . Yet, it is not necessarily the only possible solution, as the quadratic term in equation (flT))) can also be 
satisfied, (k F -[ + kp±) — ^(kp-^ + kp^kpi + k F ^) = 0. It turns out that precisely for U > U c this term has non-trivial, 
"non-balanced" solutions that are energetically more favorable than the "balanced" one and that give rise to magnetic 
order (see Appendix 1). 

Thus, equations ((9J and (fl"0|) define fcpj = k F i(U) and k F [ = k F i(U) functions which, via ([6]), give the ground 
state as a function of the external parameter U, \g) — \g(U)). This enables us to analyze the fidelity between two 
ground states \g) = \g(U)) and \g) = \g(U + SU)) in two close parameter points U and U + SU. The fidelity is then 
F(\d)(9\> \9)(d\) — an d f° r U < U c (and SU sufficiently small, i.e. SU < U c — U) we see that the fidelity is 

identical to 1 . This is a simple consequence of our mean- field approximation based on a simplified description in terms 
of single particle energy states. On the magnetic side of the phase diagram, the ground states are indeed different 
from each other, which follows from the relation k F ^(U) ^ kp\(U). In fact, from equation ^ it follows that any 
two ground states with different numbers N-\ and are orthogonal to each other. Since this is precisely the case 
in the thermodynamic limit, the fidelity between any two different ground states (in two different parameter points) 
is identically equal to zero. This is the famous Anderson orthogonality catastrophe, discussed in more detail in [f|. 
For systems with infinitely many degrees of freedom, such as those taken in the thermodynamic limit are, every two 
ground states are generally orthogonal to each other. Therefore, in order to infer the points of criticality, we are 
forced to either analyze the finite-size scaling behavi or ( see [5|), or to introduce the fidelity per lattice site and work 



directly in the thermodynamic limit (see [8j, [9( and [2l|). 

In our case though, even for finite systems the ground state is discontinuous at the point of criticality. This is due to 
the fact that the unperturbed and the symmetry-breaking perturbation Hamiltonian commute with each other, which 
results in a first-order quantum phase transition at the point of level-crossing between the ground and the first excited 
state. In other words, in the case of finite systems (TV, V < oo), the small enough changes of the parameter U — > U+SU 
will result in small changes of volumes V(k F ^) — » V(k F ^) + SV(kp^) such that |5V(A;ft)| < Vq — (2ir) 3 /V (Vb is the 
volume occupied by each one-particle state in momentum space): infinitesimal changes of the volumes V{kp-\) and 
V{kpi) are small enough to cause a change in the numbers N-\ and A^, and thus in the ground state. Therefore, for 
finite systems, the fidelity between two ground states is either one or zero - it is not a continuous function of and its 
rate of change can not be analyzed directly. 

Yet, we can use the rate of change of V(k F -[), the derivative dV< £f T ' > , to quantify the change of fidelity itself (note 
that, due to the fixed total number of electrons N, we have that dy ^ i ' > = Lengthy, but elementary 

algebra (see Appendix 1) shows that precisely at U = U c , the derivative dV fu^ diverges to infinity, thus marking 
the macroscopic distinguishability between the states from paramagnetic and magnetic phase (see Fig. [I}. 

Next, we discuss the general case of T ^ phase transitions. First, we transform the Hamiltonian HgJ to a form 
with the explicit symmetry breaking term that drives the phase transition. Using M — (Nf — N{)/2= (S z ) =X}fcWk)> 



Ground State 
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FIG. 1: (color online) The explicit dependencies fcpf = kp^(U) (black, solid line) and kpi — kpi(U) (red, dashed line) 

determining the ground state ([6]). Note that at the point of QPT, given by U c = 1, the two derivatives —Jp- and Jp become 
infinite, marking the non-analyticity in the ground state fidelity. The plot is given in rescaled quantities U — > DpU and 

kp a — > kp a /kp. 



with N a = (N a ) = J^ki^ka), n k = h k] + n kl and S% - 
^-component of the vector a of Pauli matrices, we get: 



\{&k\ - n kl ) = *t^*, where ¥ k = [ 



-fct H[ 



n SH — 



E 



UN 



■IV ) 



n k 



v 



U (N*_ 
V \~ 



M 



f2 



] and cr z is the 



(11) 



In thermal equilibrium, the state of the system is given by p — ie P(**/ H where Z = Tr[e ^^s« is 

the grand canonical partition function, /3 = X/QcbT), with fcs being the Boltzmann constant and T the absolute 
temperature and p = /u(T) is the chemical potential. Using the above expression for the Hamiltonian , we can 
write: 



eff 
SH 



k 



c. 



(12) 



Here, a k = -0E k (with E k = e k + ^ and e fe = e k - p), h z = 2/3^M and C = /3^(^ - M). Note that the 
coefficients a k , /i z and C are functions of both the coupling constant U and, through /3 and the chemical potential 
/i = jLt(T), of the temperature T as well, so that the "generalized" parameter is q = (T,U). Using the obvious 
commutation relations [n kj hy] = [hk, Sy\ — [S%, St,] = (for k ^ k'), the equilibrium state can be expressed as: 



9 Z 



,Y. k {a k n k +hJl)+C 



n* 



e a k nk e h^S k 



Tr 



n fc Tr 



(13) 



We next choose two parameter points q a = {T a , U a ) and qfo = (Tf,, Ub) defining the Hamiltonians H a — H^J (U a ) 
and Hb = Hg^J(Ub) and the corresponding equilibrium states p a = p(q a ) and pb = p(qb) respectively. The fidelity 
between the two states is then given by: 



F(p aiPb ) = ^[{p^pbp'J 2 ) 112 ] = Tr[^ 



Tr 



HaHg+fi b H b 



Z(H a )Z(H b ) 



(14) 



Using (|13[) and the expression TT[e^- aknk+h ^ s ^} = 2e Qfc (cosha/c + cosh 4^) (for the proof, see Appendix 2), the fidelity 
between two different equilibrium states p a and pb finally becomes: 



F(Pa,Pb)=Yl 



cosh a k + cosh kf- 



[cosh(a fc + ^fi) + cosh(% + [cosh(a fc - + cosh(% - 



(15) 
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with <5fc = (ak{q a ) + afc(<Zh))/2, Aak = a.k(q a ) — otk(qb), and similarly for h z and Ah z . If we choose the two points 
to be close to each other, Aak << 1 and Ah z << 1, then the fidelity can be seen as a function of otk and h z , with a 
fixed parameter difference. 

In order to evaluate the fidelity (fT5)) . we need to determine the magnetization M = M(T,U) and the chemical 
potential /i = (a(T, U), given by the pair of self consistent integral equations: 



where f(Eka) = [exp(/3Ek<r) + is the usual Fermi distribution. We used the subroutine hybrd.f from MINPACK 
[22| to solve the above system numerically. In the T — > limit, the above system reduces to equations (0 and (fTUl) 
that determine the T — ground state. The result for the magnetization is given on Fig. O The line of the phase 
transition U c — U C (T) is clearly marked, and is plotted on Fig. [3] Finally, using the numerical results for M and 
fi, we obtain the fidelity, depicted on Fig. 2] We clearly see the same line of the phase transition as the line of a 
sudden drop of F. Note that all the plots are given in rescaled quantities T — ► fc^T, U — > DpU and M — ► M/N, with 
= 0, 5t/ = 2 x 1(T 3 . We have evaluated the fidelity for ST = 2 x and SU = 0, as well as for ST = SU = 2 x 1CT 3 
and the results are qualitatively the same. 



FIG. 2: (color online) Magnetization M = M(T, U) as a function of the temperature T and the coupling constant U. The plot 
is given in rescaled quantities T — » fcsT, (7 — » DfU and M — ► M/N. 




(16) 



M = MIT, U) 



^ 0£ — 




1-0 0.0 




T 



FIG. 3: The critical line of the magnetic phase transition U c = U C (T). The plot is given in rescaled quantities T — > fcsT and 
U^DfU. 
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F-F(T.IJ) 




FIG. 4: (color online) The fidelity F = F(T,U). The plot is given in rescaled quantities T — > k B T and U — > flfU, with 
5T = 0,517 = 2 x 10" 3 . 



THE GENERAL CASE OF MUTUALLY COMMUTING HAMILTONIANS 

The above singular behavior of the fidelity indeed marks the regions of phase transitions, not only for the case of 
the Stoncr-Hubbard model but also for the broad class of systems given by a set of mutually commuting Hamiltonians 
H(q) and whose critical behavior is explained by the LGW theory, as can be seen by the following analysis. Using 
the mean Hamiltonian H = (H a + Hb)/2 and the difference AH = H a — Hf,, from equation (fl~4|) the fidelity can be 
written as: 

F(Pa, M = , . Z{H) = (17) 



Z(H + £f-)Z{H - 

Within the LGW theory, PTs occur as a consequence of the emergence of a symmetry-breaking term in the Hamil- 
tonian, given by an operator S. Thus, we can write the overall Hamiltonian as a sum of unperturbed and symmetry- 
breaking terms, H = H - hS, with h = h(q). Using [3l| Z{H) = Tr[e _/3(Ao_ ' l ' §) ] (note the implicit depen- 
dence of partition function Z on temperature T, through [3), for the first and the second derivative we obtain 
- pTr[p(H)S] = (3(S) = (3M and J^(lnZ) = X - Here, M = M(h) and x = x{h) = P((S 2 ) - (S) 2 ) are the 
generalized "magnetization" and "susceptibility" respectively. Thus, for the logarithm of the partition function (free 
energy) we obtain In Z(H)\ h + a h = h\Z(H)\ h + (3M(h)Ah + \(3x{h)Ah 2 + o(Ah 2 ). Finally, in the vicinity of the 
points of parameter space where the self-consistent field h vanishes, the fidelity (|T7|) reads as: 

F\ h = Ah ^e-^ Ah \ (18) 

According to the LGW theory, at the phase transition, the zeroth-field "susceptibility" x(0) becomes non-analytic 
and diverges. Thus it follows that the fidelity F will itself become nonanalytic and experience a sudden drop. 

When interested in the system's behavior at PTs, we usually simplify the problem by taking the unperturbed 
Hamiltonian to be constant. Yet, in general, the unperturbed Hamiltonian is also a function of the parameters, 
Hq = Ho(q). This gives the correction to the above formula for the fidelity which, introducing H — Hq — hS and 
AH = AH Q - AhS, reads as: 

F\h=Ah ^) e-toWA* (19) 

Z(H + ^)Z(S - ^ 

The "correction" term in the above product is responsible for "short-range" local correlations, while the second one 
quantifies the global "long-range" correlations giving rise to macroscopic phase distinguishability. Note though that 
even within a single phase, a system can be in different macroscopically distinguishable states - phase distinguishability 
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is not the only form of macroscopic distinguishability. Yet, it is in some sense the "extreme" version of it, which clearly 
affects the behavior of fidelity. 

Finally, we note that in the above discussion we focused on PTs driven by the local order parameter S. Therefore, 
we analyzed the Taylor expansion of F with respect to Ah deviations only, which at the second order are given by 
the generalized susceptibility x(/i). In the general case, considering the temperature deviations as well, one would 
include additional terms involving the specific heat C = C(q), again resulting in a singular behavior of the fidelity. 



BCS SUPERCONDUCTIVITY 

Next, we discuss the BCS theory for superconductivity (20j . providing us with an example of a model with mutually 
non-commuting Hamiltonians. The one-electron Bloch momentum modes are given by the fermionic anti-commuting 
operators c k(T (label a G {f , 1} represents spins with projections up and down along, say z-axis), with the one-particle 
kinetic energies taken to be, again for simplicity, e k = h 2 k 2 / (2m). The BCS superconducting Hamiltonian that 
represents the sum of one-particle kinetic and Cooper-pair interaction energies can be written in the following way 
{Vk'k = Vj* k , are the coupling constants): 

Hbcs = J2 £k£ kv dk ° + J2 Vkk ' d k>i £ -k>i d -klCkA- (20) 

ka kk' 

By hka — c ktT c k(7 we denote the one-particle number operators, while by b\ = cL.clj.1 and b k = tLfcjCfcf we define the 

Cooper-pair creation and annihilation operators respectively. Analogously to the previous case, using bk = (bk) + 5bk 
and neglecting the term quadratic in the fluctuations, we obtain the effective mean-field BCS Hamiltonian: 

H'Ms = E + n- H ) - J2(A k b{ + Alh - A* k b k ), (21) 

k k 

with A k = — Y^k> Vkk'bk' and bk = (bk)- We will use the usual assumption that the lattice-mediated pairing interaction 
is constant and non-vanishing between electrons around the Fermi level only, i.e. V kk > = — U f° r and \e k i\ < huu, 

and zero otherwise (cj.d is the Debye frequency). Using the Nambu operators [2(| T k — ip k j V>fc, where — [cjL c_fcjj 

and a is the vector of Pauli matrices (sH, the operators are given by T k = St, TZ — bk and 2T° + I = (n^ + h-kl) 
and form a su(2) algebra (see the Appendix 3). Using this notation, the Hamiltonian takes the form: 

H e Jis = E( 2e ^" - Akf k - A i f k) + Y.( £k + Akb ^- ( 22 ) 

k k 

As before, the thermal equilibrium state is given by p = L e -fmV6s~^ N ) and, using the above expression for the 
Hamiltonian H e J^, s ^ we can write: 

" PiH^s - VN) = J2 h k f k + K, (23) 

k 

where h k = {h+,h-,h° k ) = (2/?A*„ 2/3A fc , -20e k ), f k = (f+,f-,f«), K = -PJ2 k (e k + A* k b k ) and e k = e k - p(T). 

The norms of the vectors h k are given by h k = 2f3Ek, with E k = \J e\ + \Ak\ 2 ■ Similarly to what we had before, the 

coefficients h k — h k (T,V) are functions of both the coupling constant V and the temperature T, through the gap 
parameters Afc = A k {T, V) and the chemical potential \i = /i(T). Thus, we can talk of the "generalized" parameter 

q = (T, V). Since [f k , f fc /] = 0, for k f k' , we have: 

p = ± e -0(H' B "s~^ N ) = e = [lkG (24) 

z TY[e^x~ hkfk + K } Uk TT i ehkfk } 

We wish to evaluate the fidelity between two thermal states p a and pb, given for two different parameter points 
q a = (T a , V a ) and q b = (T b , V b ). Using definition (JTJ) and a k = h k (q a ) and b k = h k (qb), we have: 



Tp/~ A \ r F r/-l/2- ~l/2M/2l lr Ulfe e e e > \ /oc^ 
F[Pa,Pb) = Tr[{Pa' PbPa' ) ' \ = ~ 5 ~ • ( 2 5) 

n Jfe ('Br[e*» r *]'Iir[e I, * r *]) 1 /2 
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As for every k the operators T k form a su(2) algebra, and therefore by exponentiation define a Lie group, we can 

write e ^ LTk e bkTk e^ LTk — e 2SkTk . Also (see Appendix 3), we have that Tr[e SfcT,c ] = 2(1 + cosh^). Therefore, we 
finally have (see Appendix 3) that F(p a ,pb) — Yi k F k (p a , Pb) with [33j]: 



F k {p a ,p b ) = ^ , J , = ; V and (26) 



Tr[e SkTk ]Tr[e bkTk ] v/(l + cosh ^)(1 + cosh f) 



<oslir,=ro,l l( r£-;;x-<^l,( f'Cj \ 1 + tanh(/3 a gg)tanh(/3 b ^) £g£fc + (Afc) * J j> • «27» 



Here, we used the relation e a k =e a k - p a and a fc = 2/3°^, E% = y (el) 2 + | A£| 2 , so that cosh ^ = cosh(/3 Q ££) (and 
analogously for qt = (Vj,, ?&))■ Note the explicit dependence of all the quantities on the temperature and the coupling- 
strength, given through the superscripts a and b, denoting two parameter points q a = (T a , V a ) and = (Tf,, T4). 

Assuming that the chemical potential is also constant (p, = ef) in the region of interest, where the phase transition 
takes place, the gap parameter reduces to Afe = A, for |e| < Hlud (and zero otherwise). Thus, the self consistent 
equation for the gap Afe = — J2k' Vkk' ^~^~^^k' reads as (1 — 2f(E k >) = tanh ^ sL )'- 

f hUD tanhf^/e 2 + A 2 (T, V) , s 

1 = D F V / de 2V V =-^, 28 

J -Hun 2^£ 2 + A 2 (T, V) 

with f(E) being the Fermi distribution. 

In the T — > limit, we obtain the expression for the ground state fidelity F(\g a )(g a \, \gb)(gb\) — |(<7a|s&)|- At zero 
temperature, the chemical potential is equal to the Fermi energy ef, p(T — 0) — sf- Further, the gap equation reduces 
to A(V) — hui d / (sinh D 2 V ) — 2Hujd exp(— 2/DpV), where Dp is the density of states around the Fermi level. The 

,1/2 

zero temperature ground state fidelity is F(\q a )(q a \, |g&)(g |) = Tit -4s I 1 H (Ek-er) +A(y a )A(v b ) 

* 6 J \\yai\ya\, iyo/\ya\j llk^y v / (e fc -£F) 2 +A(y a ) V(efe-E-F') 2 +A(V 6 ) 2 

which matches the expression one obtains for T = 0. In other words, we see that, as in the case of the Stoner- 
Hubbard model, the point of criticality of the T — QPT can be inferred from the mixed state fidelity between the 
thermal states. The phase transition from superconductor to normal metal happens at V — 0. Thus, for V a — > 0+ 
and Vb = V a + SV — > SV > 0, the fidelity between the corresponding ground states of Fermi sea \gp) and the BCS 

( _ i \ 1/2 

superconductor \gscs) becomes: \(gF\gBCs}\ = TL -4s i + / le "~ e ■ We sec that the BCS model at 

V-! y ^/(e fc -e F ) 2 +A(i5y) 2 J 

T = features the Anderson orthogonality catastrophe, just as the Stoner-Hubbard model. 

As in the previous case, in obtaining the numerical results for the fidelity, we used the subroutine hybrd.f from 
MINPACK [22|. Again, all the numerical results are given in rescaled quantities T — > fc^T '/ '(hwo), V — > DpV and 
A — > A/(fkJr>). The result for the gap is given on Fig. \5\ The line of the phase transition is clearly marked as the line 
along which the gap becomes non-trivial, and is presented on Fig. [S] Finally, the fidelity, with 8T = 0,5V — 10 -3 , 
is plotted on Fig. [71 We varied the parameter differences ST and SV, and all the results obtained show the same 
qualitative picture the fidelity exhibits a sudden drop from F = 1 precisely along the line of the phase transition. 

As already discussed, in the case of mutually commuting Hamiltonians, the fidelity reduces to the quantity 

C(p ai Pb) = , . Z{H) . (29) 
Z(H + Af)Z(H - Af) 

which, through relation (|18[) . establishes the connection between the singular behavior of the fidelity and the corre- 
sponding susceptibility (or the heat capacity, etc.). In the non-commuting case, the same relation between C(p a ,pb) 
and x is still valid (see Appendix 4), yet the fidelity is in general not identically equal to C. 

If we approach a line of the phase transition along a curve q = q(a), a € 1, we find: F — C(l + J-da 2 ), where 
da defines the difference Sq. For reasons of simplicity, we omit here the explicit, lengthy expression for T . Although 
relatively complicated and difficult to study directly, it is evident that, apart from the Hamiltonian's eigenvalues, it 
is also explicitly given by the rate of change, with respect to a, of the Hamiltonian's eigenbasis. This is also evident 
from the fact that for the case of mutually commuting Hamiltonians, when the eigenbasis is common, C = F and 
thus T = 0. In the commuting case, the change of state is given by the change of the Hamiltonian's eigenvalues 
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A =A(t, V) 




FIG. 5: (color online) The gap A = A(T, V) as a function of the temperature T and the coupling constant V . The plot is given 
in rescaled quantities T — > fesT /(fk/jo), V — > DfV and A — > A/(hcuD)- 



FIG. 6: (color online) The critical line V c = V C (T) as a function of the temperature T. The plot is given in rescaled quantities 
T -> k B T/{kuj D ) and V" -> DfV . 



only, while in the non-commuting case, the overall change of state is given by the change of both the eigenvalues 
and the eigenvectors. Thus, the singularity of T can on its own mark the drastic change in the structure of the 
system's eigenbasis, thus bringing about the finite difference between C and F. This is indeed the case for BCS 
superconductors, where the difference C — F becomes non-trivial precisely along the line of the phase transition, 
where C < F, see Fig. [3] We see that, as intuitively expected, the state of a system exhibits an abrupt change in its 
structure along the line of the phase transition both in terms of its eigenvalues, as well as in terms of its eigenstates 
(for the structural analysis of the system's eigenstates given by a parametrized Hamiltonian H(q), see for example 
[23| and Q). 

Another way to quantify the structural change of the eigenvectors is through the Uhlmann connection and the 
corresponding mixed-state geometric phase (l9| . the mixed-state generalization of the Berry connection and phase. 
The fidelity F(p a ,p b ) = Tr-y/ \Jp~ a pb\Jp~a can also be expressed as F(p a ,pb) = Tr [IVPoVpfel] i where \A\ = (AA^) 1 / 2 
represents the modulus of the operator A (see [IBj]). The Uhlmann parallel transport condition (i.e. the connection) 

Tr \fp\\fpbUab (see equation (9) from [25]). In other 



is given by the unitary operator U a b, such that F(p a ,pb) 

words, the connection operator defining the parallel transport is the inverse of the unitary V a b, given by the polar 
decomposition _lj] \f~p\\f~pb — \\/ r Pa\ /r Pb\Vab, i.e. U a b = V^b- Note that the parallel transport condition (the connection), 
given by U a b, induces both the local Uhlmann curvature two- form, as well as the global mixed-state geometric phase 



F-F(T,V) 
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FIG. 7: (color online) The fidelity F = F(T, V) as a function of the temperature T and the coupling constant V. The plot is 
given in rescaled quantities T -> k B T/(hujD) and V -> DfV, with <5T = 0, 8V = 10~ 3 . 



C(T, V)-F(T, V) 



FIG. 8: (color online) The difference C(T, V) — F(T, V) as a function of the temperature T and the coupling constant V . The 
plot is given in rescaled quantities T — » fcsT/ (HcJn) and V — » DfV, with 5T = 0, SV = 10~ 3 . 



(see equations (11) and (12) from [25[). Let us define 



H(p ai p b ) = Tr 



(30) 



Obviously, in the case of mutually commuting Hamiltonians if = i 1 (= C) and the Uhlmann connection is trivial, 
Uab = I- In the case of the BCS model, we have evaluated the quantity if — F, and the result is qualitatively 
the same as for C — F. In other words, along the line of the phase transition we have found the strict inequality 
if — F = Tr \y/p\y/pb\(U a b — i) < const < 0. From the Cauchy-Schwarz inequality it follows then that U a b — 1^0, 

or U a i, ^ I — a clearly abrupt change of the connection operator U a f, occours and it becomes non-trivial in the vicinity 
of the line of the phase transition. Since the two parameter points q a and qi, are taken to be close to each other, such a 
behavior implies the non-analyticity of the local Uhlmann curvature form along the line of the phase transition, which 
in turn results in generally non-trivial global Uhlmann mixed-state geometric phase (see for example the discussions in 
[l3l ] and [3])- This can be seen as a mixed thermal state generalization of recent results [lj| on the relation between 
QPTs and Berry geometric phases. Further, we have that along the line of the phase transition we have C < H < F, 
while C ~ if ~ F ~ 1 otherwise. 
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CONCLUSION 



In this paper we analyzed the fidelity approach to both zero temperature (quantum) as well as finite temperature 
phase transitions. It is based on the notion of quantum state distinguishability applied to the case of macroscopic 
many-body systems whose states determine the global order of the system and its phase. We focused on the two 
particular cases of the Stoncr-Hubbard model for itinerant electron magnetism and the BCS theory of superconduc- 
tivity, as representatives of two distinct classes of physical systems: those defined by a set of mutually commuting 
Hamiltonians and those defined by mutually non-commuting Hamiltonians, with respect to a given parameter space. 
We found that in both cases the fidelity can mark the regions of PTs by its sudden drop in the vicinity of the transition 
line. We discussed in detail the general case of mutually commuting Hamiltonians, where in the case of symmetry- 
breaking LGW type of transitions the non-analyticity of fidelity is a direct consequence of the non-analyticity of the 
corresponding "generalized" susceptibility. The case of mutually non-commuting Hamiltonians is more complex as 
there the structure of the Hamiltonian's eigenvectors directly affects the system's thermal state thus introducing a 
novel feature relevant for the macroscopic phase distinguishability. On the example of the BCS superconductivity we 
showed that the phase transition is accompanied by the abrupt change of the Hamiltonian's eigenvectors, quantified 
by the differences C — F and H — F. The second one is of particular interest as it determines the Uhlmann connection 
and the associated mixed-state geometric phase, which generically becomes non-trivial in the vicinity of the phase 
transition only - a mixed state generalization of the emergence of non-trivial Berry geometric phase in the vicinity of 
criticality for QPTs. 

There are two main extensions of this work. First, the study of more general phase transitions that fall outside the 
standard symmetry-breaking LGW paradigm. Second, the general study of the structural analysis of the system's 
eigenstates in the case of mutually non-commuting Hamiltonians, given within the framework of the mixed state 
fidelity and the Uhlmann geometric phase. 
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APPENDIX 1 



In this Appendix we prove various technical statements relevant for the analysis of the ground state fidelity in the 
Stoner-Hubbard model. 

• Proving the relation ef = £fe FT = £k Fl - 

Using N a = J^" dN ai the ground state energy can be written as: 

E 9 = J o E k] dN ]+ J^ E kl dN l -U-^= J Ek dN^+J Ekd^+U-^. (31) 

Since it is the local minima, for given N — JVf + Aj_ we have that: 

d(E g -»N) N i 

— — — = £fcFT _ M + C /_ = , 

d{E g -nN) ll + U N ^-0 

— _ £kFl -^ + U —-0, (32) 

and we immediately get £fc FT = £k F ] + Un^ = £k F [ + Un^ = Ek Fl = £f, with fi — ef, for T = 0. 

• Deriving U c from Stoner criterion: 

The density of states is defined as D{e) — y^JT m ^ ne paramagnetic phase, where e is the one-particle energy 
and N„{e) is the number of states for each spin direction whose energy is smaller or equal than e. For each 
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spin projection, we have N a /V — k F /((>Tr 2 ), while from e(fc) = j^k 2 we have that k = 



NAe) = 1 

V 6TT 2 



(2m e )3/2 and . 



D(s) = 



1 dN„(e) 
V de 



3 1 
26^ 



2m 



3/2 



r _ 3 N„ _ 3 N 



Thus, 



(33) 



The difference between the up and down single particle energies and E k [ defines the self consistent magnetic 



field A 



2U 



M (A) 



The phase transition occurs when A becomes non-zero, leading to the equation 1 = 



~ " "" — |a=0j which also gives the vanishing of the susceptibility denominator in the random phase approximation 



V dA I 

(Stoner enhancement factor) [20j. One has then U c Dp 



1, and therefore we find U c = 



4_ V 



• Deriving U c from ([9]) and (fTOj): 

Using x = y — kj?i, equation 



can be rewritten as f(x, y) = a; 3 + y 3 — a = 0, with a = 6tt 



2N _ 



Also, the quadratic term from equation (|10|) that is for U > U c responsible for non-trivial magnetic solutions 



2/r _ 3 



k F Ur 



reads as g(x,y; U) = (x + y) — ^(x + xy + y 2 ) = 0, with a — 3tt 2 ^- 

First, we show that for U > U c (see equation ([8|). curves / and g have intersection in two points, symmetric 
with respect to the y = x line, while for U — U c they touch precisely in the point (x,y) — (k F ,k F ). In the 
coordinate system (x,y), rotated by the angle ip = ir/4 and translated by the vector (A,0) from the system 
(x,y), the curve g(x, y; U) = becomes [3J|: 



A 2 B 2 



= 1, 



(34) 



with A = 



and B = yf^ 1 - a real ellipse with the main axes A = A(U) and B — B(U). As the 

parameter U increases from (when A,B — > oo), the main axes A and i3 decrease and eventually the ellipse 
g touches the curve / in the point (x,y) = (2 A, 0), which in the original coordinate system correspond to the 
point (x,y) — (kp,kp). In other words, 2A — \/2kp and U — U c . Further increase of U beyond U c results in 
even smaller ellipses that intersect the curve / in two symmetric points (note that the curve / is itself symmetric 
along the y — x line) . 

To show that for U > U c the "balanced" solutions of © and (|10[) are indeed the physical ones, we calculate the 
total ground state energy E g of the system: 



E g = V f F1 £ k] D{E k] )dE k] +V f ^ E H D(E H )dE kl -VUn ini . 

J -Bot J Eqi 

Upon evaluating the above integrals, using M = (JVj — N{)/2, the total ground state energy reads: 



v h 2 

20tt 2 m 



9 N 

37r V 



5/3 



2M\ 5/3 / 2M\ 



5/3' 



UN 2 
AV 



(35) 



(36) 



OE 

From the above expression, we see that for every value of the coupling constant U, -arf(U, M — 0) = 0, while 



d A E a 



gjtf %~(U, M — 0) changes the sign from positive to negative precisely in U = U c . Thus, we indeed have that for 
U > U c the "balanced" magnetic solutions with M ^ have lower energy. 

• Deriving the maximum fidelity, i.e. the maximum derivative dkp/dU for U — U c , from ([9]) and 

god: 

From df = and dg = 0, we get (|| = f x , % = x, etc.): 

fxx + fyij = 0, 

gxx + g y y = -g v . (37) 



Solving it for x and y, we obtain: 



x 

y 



9-u 



fxgy fygx 



fy 

~~ fx 



(38) 



14 



Finally: 

'Y- 1 3 h £ T-l-7/1 T <j;2 

(39) 



3 kf x + y 1 
A DpV 2 x — y xy 



y 2 

2 

— X 



Thus, in the limit U — > C/ c + (note that g(x, y; U) — is valid for U > C/ c ), we have that x — » Z/ — > &f— and 
i — > +oo, y — * — oo (or vice versa, for the x < y solution branch). In other words, ^—%p^ — 47r(fcf^) 2 = 
4irx 2 x —>■ +oo and dV ^^ — * — oo, or vice versa. 

APPENDIX 2 

Below, the basic identities of the su(2) algebra of the electron spin operators leading to equation (|15|) are proven. 
We discuss operators having a fixed momentum and drop the index k for convenience. 

• S'-algebra identities: 

In general, operators defined by S a = & a S^pCp, where cj,, and cp are either bosonic or fermionic operators and 
S^p is a matrix representation of a spin 5 algebra, i.e. satisfying [S a , S b ] a p — iE ahc S°p, define a su(2) algebra. 

This is the case of the electron spin operators defined by S = where ^ = [cj cj] and a is the vector of 

Pauli matrices, and given by S z — |(rif — hi), S + = cjcj, and 5~ = cjcf (with S' ± = S" 21 ± iS' 3 '). 

Using the anti-commutation relations for the one-electron modes cj., one can easily obtain the following com- 
mutation and anti-commutation identities: 

[S*,^] = ±S ± , [S + ,S-] = 2S Z , (40) 

{s z ,s ± } = o, {s + ,s-} = i s . 

with I s = n| + h[ — 2n^hi = (2 — n)n, n = + rjj,, the projector onto the subspace of the single-occupied 
states, being an identity operator for this algebra, since I 2 = I s , I S S = SI S = S and (S) 2 = |/ s . One of the 
immediate consequences of the above commutation relations is that {S x , S y , S z } form a su(2) algebra. 

Further, we can obtain the traces: Tr[/] = 4, Tr[I s ] = 2 and Tr[5] = 0. 



• Evaluating e s and Tr[e 



+oo 1 ^ +oo .. 

n! 1 ' * — ' n 

n— n— 1 



£^) n ='+E^)" ( 41 ) 



For (hS) i , we have: 

(K<?) 2 = [i(/i+S- +/ l -S' + ) + /i z S' z ] 2 (42) 
= J(^ + ) 2 (5-) 2 + i(^) 2 (5+) 2 + h 2 z (S*) 2 + \(h+h-){S~,S+} + \{h z h+){S*,S-} + \(h z hr){S\S+}. 

Using the 5-algebra identities (|4T)]) . we get: 

(h"sf = \{Kfi s + \{h+h-)i s = \\h\ 2 L (43) 

where (3?) |/i| 2 = (/i z ) 2 + = h 2 . Thus, we obtain: 

+oo /i \ 2n +oo / , \ 2n 



15 



{I -Is 



+00 J 
n=0 y ' v 



+00 ^ 

^ (2n + Y)\ \ 2 

.n=0 v ' v 



2n+l 



(/-/.i+<osh('^I s + 2 S inh(^M. 



(44) 



Using the trace formulas, we finally obtain: 



Tr^" 6 ] =2 1 + cosh - 



(45) 



Evaluating e afl and Tr[e Q "]: 

Enlarging the su(2) algebra of the S operators to the operator h = h^ + hi it is useful to know its algebraic 
properties and their relation to the other operators. We first evaluate h k . Since h\ = h a and h 2 = (h^ + n|) 2 = 
h + 2n^ni we find that h 3 = —2h + 3h 2 and in general h k = —dkh + bkh 2 , where at and bk satisfy the recurrence 
equations at+i — 2bk and bk+i = 3bk — ak, leading to: h k = — (2 fc_1 — 2)h + (2 fe_1 — l)h 2 . The exponential e an 
is given by: 



I 



,2ct 

2e a + - I n 

2 2 



-,2q 



/ + un + ufi 



(46) 



which can be easily verified since h = + can only take the values n = 0, 1, 2. 

Using the obvious trace relations Tr[n] = 4 and Tr[n 2 ] = 6, we also get Tr[e Qn ] — 4 + Au + 6w, or: 

Tr[e Qft ] = (e Q + l) 2 . 



(47) 



The number operator commutes with the S operators, and one has nl s — I s n — I s and nS — Sn = S. Finally, one 
obtains: 



Tr 



^{an+hS) 



l + e c 



2 e a ( cosh a + cosh ■ 



(48) 



APPENDIX 3 

In the following, we obtain the basic formulas for the su(2) algebra of the Nambu operators T used in evaluating 
equation PB"]) for the fidelity between two states of a BCS superconductor in thermodynamical equilibrium. 

• T-algebra identities: 

We now consider the operators defined by T = ip^'ip, where $ — [cj cjj and a is the vector of Pauli matrices, 

and given by f z = |(n T +nj. - 1), f+ = c\c\ and f~ = c x c T (with f ± = f x ± if«). These operators f a are of 
the form considered before, since the fermion anti-commutation relations are invariant under the electron-hole 

transformation, i.e. under the interchange c «-> & . Similarly to the previous case of the S-algebra, we obtain 
the commutation and anti-commutation relations: 

[f °, f±] = ±f ± , [f + ,f-\ = 2T°, (49) 

{f°,f ± } = o, {f+,f-} = i t . 

with I t = 2hfhi — (fij + hi) + 1 = (n — l) 2 , the projector onto the subspace of the empty and doubly occupied 
states, being an identity operator for this algebra, since If — I t , I t T = TI t = T and (T) 2 = |j t . One of the 
immediate consequences of the above commutation relations is that {T°, T + , T~} form a su(2) algebra. Further, 
simple algebra gives: 



Tr[f] = 0, Tr[/ t ] = 2. 



(50) 
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Note that in the above definition of $ , the operators cj and cj, have the opposite momenta, k and —k, respec- 
tively. Had the momenta been defined to be the same, the two sets of operators S and T would commute with 
each other, i.e. [S a ,f b ] = 0. Also, one would have I s + I t = I. 

Evaluating e sf and Tr[e af ]: 



+00 1 _ +00 



(51) 



Similarly to the previous case, for (of ) 2 , we have: 



(aff = [a z f° + l -(a+f- + a~T +)f = \\a\ 2 I u 



(52) 



where \a\ 2 = (a z ) 2 + a + a = a 2 . Thus, we obtain (Z t = I — It 



+00 



= (/-/t) + 



£(2n+l)l 



1 /a\ 2 ™ 
(2^)! 15 

n— v J 



y 1 

^ (2n+ 1)! V2/ 



2n+l 



-(of), 
a 



e^Z t+ cosh^)i f+ 2sinh(|) M. 



(53) 



Using the trace formulas (|50|) . we finally obtain: 



Tr[e^] = 2(l + cosh^) 



(54) 



Evaluating e 2 ^ = e and Tr[e CT ]: 

From the above formula (fB"4"]) and the known trigonometric relation (cosh^-) 2 = i(l + cosh a), we have that 



Trfe 57 ] = 2 (l + cosh |) = 2 ( 1 + y i(l + cosh c)J . Thus, by evaluating e 25T we obtain the result for coshc 

and therefore Trfe 57 ]. First, we obtain the general expression for e aT e bT . Using the result (fS"3"|) for e aT (and 
analogously for e bT ), we get: 



e ST e bT =Z t + [ cosh - cosh - I I t +2 



cosh — sinh — I - + ( sinh — cosh - — 
2 2) b V 2 2a 



T+4 sinh -sinh- .(55) 



In deriving the above expression, we have used the above identities (|50|) . As in (|52|) . using (|49|) we obtain: 



(ST)(bT) = - y{ab)I t + 2i{ax b)T 



(56) 



Therefore, we have: 



e ST e bT 



( u a , b \ ( ■ 1 « • 1 h \ (« & ) 

cosh — cosh — + sinh — smh - - 

\ 2 2 J V 2 2/06 



. , a , b\ a ( a . b\ b (a x b) 

smh — cosh - — M cosh — smh - — \- 1 

2 2 i a I 2 2 b ab 



(57) 



Applying the above result twice, we finally obtain the expression for e 



2cT. 



e 2cT = e f T e& T e f T = ^ 



cosh ^ cosh ^ 
2 2 



sinh — sinh — | ^—1. 
2 2/ a& 



sinh — cosh — | — + 2 ( sinh — | - + ( cosh 1^ ( sinh , 

2 2/ a V 2) b \ 2 / V 2 J a 2 b 



b\ (ab)a 



T. 
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(58) 



Comparing the above result with the expression (|53p , we eventually end up with the expression for cosh c: 



cosh c = ( cosh — cosh — 
2 2 



sinh — sinh — | — — 
2 2 ab 



(59) 



which reduces to (|27|) for particular values of = hk(q a ) and 6fc = /ife(g , l 



APPENDIX 4 



In this Appendix we prove that in the case of mutually non-commuting Hamiltonians a relation analogous to (jTSJ) 
holds between C, given by equation (|29j) . and the susceptibility x- As in the commuting case, for simplicity, we 
consider a Hamiltonian H = Hq — hS with the symmetry-breaking term S, and h — h(q). Note that in this case, the 
two terms in the Hamiltonian do not commute with each other, [Hq,S] =/= 0. Thus, we have the following imaginary 
time Dyson expansion around the point h = [36| : 



-P(H -hS) 



e-P"°+h f dTe- p "°S{T)+h 2 f dr [ d Tl e~^« S(t)S( Ti ) \ , 
Jo JO JO 



(60) 



with S(t) = e TH °Se T °, From the above equation, we obtain the expressions for the magnetization M — (S) and 
the susceptibility \ = given by derivatives of the partition function Z. First, the magnetization can be expressed 
as (using the commutativity between the partial derivative and the trace, ^-Tr[-] = Tr^[-]): 



1 d\nZ 1 1 dZ lid 



1 1 



13 



1 







,-/3H 



M = ^ = J-zm = T-zHh^-^ = Ji L ^-°"°^-"-] = - f I dr-nl—S] = (S>.(6i, 

The susceptibility is then: 



X 



8M d ( 1 1 dZ 



1 1 d 2 Z 11 (dZ 



dh dh\l3ZdhJ (3 Z dh 2 (3 Z 2 \dh 
The second term is obviously equal to (3M 2 — (3(S) 2 , while the first term can be transformed as follows 37J: 



(62) 



1 1 (d 2 Z 



1 1 







1 



de-? H 
dh 



S} = - dTTr[e-P 6o S(T)S]= (Lt(S(t)S). (63) 



Thus, the susceptibility is given by x = f£dT[{S(T)S) - (S) 2 }. From this, the Taylor expansion for Z, 



Z Z n \l + (3Mh + ^[3 2 M 2 h 2 + ^[3 X h 2 



(64) 



is identical to the one obtained for the case of mutually commuting Hamiltonians, and therefore a relation analog to 
equation (fTS|) holds between C and the susceptibility x- 
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Fidelity is closely related to various distance measures on the set of quantum states: the closer two states are, the more 
similar (less distinguishable) they are (for an overview of different distance measures induced by the fidelity, see [3]). 
The use of the fidelity approach to thermal phase transitions was also suggested in 
For the proof of the relation ef = Ek F1 = Ek Fi see Appendix 1. 

In the following, we discuss the case of a 3-dimensional spatial lattice. This treatment could be easily generalized to a 
different lattice dimensionality. 

Strictly speaking, for U > U c we speak of multi (two- fold) valued functions kpf(U) and kpi(U). 

For reasons of simplicity, we use formulas for canonical rather than grand canonical ensemble. Introducing chemical 
potential results only in a "shift" of the one-particle energies from the unperturbed Hamiltonian Ho- 
Note the similarity, and the difference, to the previously introduced operators ^fc. 



We use the well known identity cosh ^ = y i(l + coshcfe). 

Alternatively, we could first translate the coordinate system by the vector {\/2A, V2A) and then rotate it by the angle 
<p = 7r/4 to obtain the same system (x,y). 

Using A = A x ± iA y and A = A z , we can write A — A x e x + A y e y + A z e z — A + e+ + A~e- + A°eo, where e± = e x ± ie y 
and e*o = e z (and analogously for B). The scalar product between two vectors is then A ■ B = A X B X + A y B y + A Z B Z = 
\{A + B~ + A~B + ) + A°B°. Therefore, A\ = A + A~ + (A ) 2 . 

We consider the particular case of the h — expansion, but all the results obtained can be easily generalized for the h ^ 
case. 

In order to avoid a possible confusion, here we explicitly denote that the derivatives are taken for h = 0. 



